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ABSTRACT 



Context, X-ray bright active galactic nuclei represent a unique astrophysical laboratory for studying accretion physics around super- 
massive black holes. 

Aims. 4U 1344-60 is a bright Seyfert galaxy which revealed relativistic reflection features in the archival XMM-Newton observation. 
Methods. We present the spectroscopic results of new data obtained with the Suzaku satellite and compare them with the previous 
XMM-Newton observation. 

Results. The X-ray continuum of 4U 1344-60 can be well described by a power-law component with the photon index » 1.7 modified 
by a fully and a partially covering local absorbers. We measured a substantial decrease of the fraction of the partially absorbed 
radiation from around 45% in the XMM-Newton observation to less than 10% in the Suzaku observation while the power-law slope 
remains constant within uncertainties. The iron line in the Suzaku spectrum is relatively narrow, cr = (0.08 ± 0.02) keV, without any 
suggestion for relativistic broadening. Regarding this, we interpret the iron line in the archival XMM-Newton spectrum as a narrow 
line of the same width plus an additional red-shifted emission around 6. 1 keV. 

Conclusions. No evidence of the relativistic reflection is present in the Suzaku spectra. The detected red-shifted iron line during 
the XMM-Newton observation could be a temporary feature either due to locally enhanced emission or decreased ionisation in the 
innermost accretion flow. 

Key words. Galaxies: active - Galaxies: Seyfert - Galaxies: individual: 4U 1344-60 
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1. Introduction 

Bright active galaxies provide a good opportunity to probe 
accretion physics in the strong gravity regime and to reveal 
the basic parameters of the accreting black hole (see e.g. 
lRevnolds~ & Nowak l2003h . In particular, the black hole angu- 
lar momentum can be constrained from the relativistic fea- 
tures in the spectra originating in the closest neighbourhood 
to the black hole. Suitable tools for this purpose are the high- 
throughpu t X-ray detectors on space mis sions such as XMM- 
Newton dJansen et alj I20Q lb and Suzaku dMitsuda et alj I2007) 
whose data analysis we are presenting in this paper. 

A sizeable sample of active galactic nuclei (AGNs) 
with relativistic ally b roadened iron lines was studied by 
iGuainazzi et alJ {2006) and later by Ide La Calle Perez et alJ 
d2010t) . Different samples of Seyfert galaxies ob served by the 
XMM- Newton satellite were also exa mine d by Nandra et al 



, iBrenneman & Reynolds! d2009t) . and Bhavani & Nandra 
Broad iron lines in quasars were studied by 



Jimenez-Bailon et al. (2005). All groups concluded that X-ray 
spectra of a substantial fraction of AGNs possess a relativisti- 
cally broadened iron line. The fraction increases with longer- 
exposed observations providing better statistics for spectro- 
scopic studies. 

4U 1344-60, the subject of the study discussed in this pa- 
per, is a bright (about 2 millicrabs) and nearby source but it 
has not been intensively studied so far due to its low Galactic 
latitude (|Z?|=1.5°). Although its X-ray spectrum is heavily ab- 



sorbed by the Galactic interstellar matter (Nh ~ 1 X 10 cm" ) 
it has clearly revealed an emission excess in the broadened red 
wing of the iron line in a previous XMM-Newton observation 
dPiconcelli et al.ll2006l) . 

4U 1344-60 was discovered for the first time by the 
X-ray satellite Uhuru. In the XMM-Newton observation of 
Centaurus B, 4U 1 344-6 appeared at the e dge of the 
field (Obs.ID: 00921 40 lOlV lPiconcelli et al] d2006l) studied the 
EPIC-pn spectrum together with the optical observations of the 
source. The optical spectrum allowed for a classification of the 
source as an active galaxy of an intermediate type. The cosmo- 
logical redshift was derived to be z = 0.012 + 0.002. Its flux, 
fi-AOkeV ~ 3-6 x 10" 11 erg cm" 2 s _1 , makes it one of the bright- 
est AGN in the hard X-ray sky. The luminosity, £f^[$k eV = 
1.5 x 10 43 erg/s, is ty pical for Seyfert g alaxies. 4U 1344-60 was 

aU i2010l) 



observed by Integral dBird et 

/2notv=( 3 - i8±o -° 8 )>< io " i! 

(4.1 



with its hard X-ray flux 



erg cm 2 s 1 



and f Inte g ml 
a " u ■/ 40-100 keV 



0.2) x 10 n ergcm 2 s 1 . It was also detected in the Swift- 



(9.0 + 0.3) x 10" 11 erg cm" 2 s" 1 
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BAT survey with ff^Cv 
(ICusumano et al.ll2010h . 

Residuals in the energy range, where emission fea- 
tures by iron l i ne are expected, were formerly analysed by 
iPicon celli et alj (120061) in the XMM-Newton observation. A 
model consisting of the power-law continuum plus a narrow un- 
resolved Gaussian line gave an unacceptable fit. If the line profile 
was allowed to be distorted by relativistic effects the fit was sig- 
nificantly improved. The inner disc radius was found to be within 
10 gravitational radii (r g = The outer one was larger than 
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Fig. 1. XISO, XIS 1, XIS3, and HXD/PIN (all background sub- 
tracted) light curves in the 0.5 - lOkeV, or 15 - 60keV energy 
range, respectively. The time bin size is 4096 s. 

90 r g with a fixed radial emissivity index q = 2.5. The inclination 
of the disk was found to be around 20 degrees. The equivalent 
width of the line was about 400 eV. 

4U 1344- 60 also exhibited a com plex absorption in the X- 
ray spectrum. iPiconcelli et al.l d2006l) suggested the presence of 
a partially covering absorber to explain a remarkably flat power 
law describing the data. A combined analysis of the XMM- 
Newto n and Integral observation was later done by Panessa et al. 
(2008) who concluded that the spectrum is affected by one fully 
and two partially covering absorbers with substantial column 

densities (A^ cl * 5 x 10 22 cirT 2 , and N^ 2 ~ 4 x 10 23 cirT 2 ) 
and covering fractions around 50%. 

In this paper, we present new spectroscopic results of the re- 
cent observation of 4U 1344-60 obtained with the Suzaku satel- 
lite. We compare them with the archival XMM-Newton data, 
which we re-analysed. The paper is organised as follows. Sect. [2] 
describes the data reduction of the Suzaku XIS and HXD/PIN, 
and also the re-analysis of the XMM-Newton data. The ba- 
sic timing and spectral properties of the Suzaku observation 
are presented in Sect. [3] We compare the spectral results from 
the Suzaku observation with the archival XMM-Newton data in 
Sect. |U The achieved results are discussed in Sect. [5] and the 
main conclusions are drawn in Sect. [6] The timing and spectral 
analysis of Centaurus B is presented in the Appendix A, as well 
as the estimate of its contribution to the 4U 1344-60 HXD/PIN 
spectrum. 

2. Data reduction and analysis 

4U 1344-60 was observed by the Suzaku satellite on January 
201 1 (obs. ID 705058010) for the total exposure of 100 ks. The 
observation was performed using X/5-nominal position. It was 
accompanied by a short lOks observation pointed to the radio- 
galaxy Centaurus B (obs. ID 705059010), which is a neighbour- 
ing source on the X-ray sky. Due to its proximity, the HXD/PIN 
spectrum of 4U 1344-60 is likely to be contaminated. The short 
observation of Centaurus B was performed to measure its flux 
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Fig. 2. XISO (background subtracted) light curves in different 
spectral bands: 0.5-2keV (upper panel), and2-10keV (middle 
panel). Their ratio is plotted in the bottom panel. The time bin 
size is 4096 s. 



and spectral properties below lOkeV with the XIS detectors, in 
order to estimate its contribution to the HXD/PIN spectrum of 
4U 1344-60. This observation strategy was possible only thanks 
to the very low short-time variability of the Centaurus B and 
its simple power-law shaped X-ray spectrum, as was suggested 
from the previous observations by the ASCA and XMM-Newton 
satellites (see Appendix). 

We used the Heasoft package version 6. 1 1 . lQ for the data re- 
duction and also for the subsequent spectral and timing analysis. 
The data were processed standardly following the Suzaku Data 
Reduction Guided (version 4). For all XIS detectors, we com- 
bined both 3x3 and 5x5 modes to extract the event files. The 
source spectra of 4U 1344-60 and Centaurus B were obtained 
from a circle around the centre of the point spread function with 
the radius of 260 arcsec. We defined the background extraction 
region as an annulus around the source circle with the outer ra- 
dius of 360 arcsec to avoid any contamination from the calibra- 
tion source. We created the related response matrices and ancil- 
lary response files using the tools xisrmfgen and xissimarfgen. 
The HXD/PIN spectra were reduced with the tool hxdpinxbpi. 
The tuned background files were used for modelling the non 
X-ray background. The standard method was used to account 
for cosmic X-ray background as described in the Suzaku Data 
Reduction Guide. 

Cross-normalisation between the Suzaku detectors were 
fixed to 1 for XIS 0, free for XIS 1 and XIS 3, and fixed to 1 . 16 for 
HXD/PIN. Fits were performed in the l-10keV and 15-60keV 
energy ranges for the XIS and HXD/PIN, respectively. The data 
below 1 keV were not included due to large Galactic absorption. 
The XIS 1 data in the energy range 1 .6-2.3 keV were ignored due 
to calibration uncertainties. We used C-statistics (ICashl [19791) 
for fitting the (unbinned) data, which employs the appropriate 
Poisson distribution of count fluctuations. However, we express 
the goodness of the fit also with the more familiar^ 2 statistics. 
Only for this purpose, we binned the spectra to contain at least 
30 counts per bin because the Gaussian smooth distribution of 
count fluctuations is good only for a sufficient number of counts 



http://heasarc.nasa.gov/lheasoft/ 
http://heasarc.nasa.gov/docs/suzaku/analysis/abc/ 
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Fig. 3. XISQ (black), XISl (red), XIS3 (green), and HXD/PIN 
(magenta) spectrum. The underlying model is an absorbed power 
law with Fx 1.7. The spectrum is binned for plotting purposes 
only. 
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Fig. 4. Contour plot of the rest energy of the iron Ka line and 
its width as measured with Suzaku. The individual lines corre- 
spond to 68, 90, and 99% confidence levels. The best-fit values 
are indicated by a small cross inside the contours. 



per bin (see e.g. lNousek & Shudl989L for comparison of C- and 
X 2 - statistics). The quoted errors in the text correspond to a 90% 
confidence level for one interesting parameter. 

Because our spectral re sults with the Suzaku o bservation dif- 
fer from that reported by iPiconcelli et alj d2006l) we also re- 
analysed the archival XMM-Newton data and performed a si- 
multaneous fit to both data sets. We used the recent version 
of the SAS software (ll.O.Ofl 4U 1344-60 was observed with 
XMM-Newton in August 2001 (obs. ID 0092140101) for the to- 
tal exposure of 37 ks. The observation was aimed to Centaurus B. 
4U 1 344-60 is well detected in the EPIC-pn field-of-view while 
it is on the edge of the EPIC-mosl and outside of the EPIC- 
mos2 field-of-view due to a relatively large offset 14 arcmin). 
Therefore, we used only the EPIC-pn spectra. 

The data were standardly screened - only patterns 0-4 were 
used, high background flaring events were ignored (background 
rate > 0.4 ctss~'). The source was defined as an ellipse (in 
the sky xy-coordinates: 18800.5, 12880.5, 1946.4, 909.7, 327), 
which better corresponds to the shape of the p oint spread func- 
tion due to its offset position (Read et al. 201 1). The background 
extraction region was defined from a nearby circle in the source- 
free region of the chip. The results of our re-analysis of the 
XMM -Newton EPIC-pn s pectrum are consistent with the results 
bv lPiconcelli et alj d2006l) . 

3. Suzaku view of 4U 1344-60 

The time-averaged flux of 4U 1344-60 observed by Suzaku is 
similar to the flux measured by the XMM-Newton satellite, i.e. 
/2-iokeV ~ 4 x 10~ n ergcirT 2 s _1 , or the luminosity L2-i0keV ~ 
1.2xl0 43 erg s _1 , respectively. During the observation the bright- 
ness of the source increased by a factor of 2 - see Figure[T]where 
background subtracted light curves of all Suzaku detectors are 
shown. However, the spectral shape did not change significantly, 
as it is suggested by the almost constant hardness ratio defined 
as a ratio of 0.5 - 2keV to 2 - lOkeV spectra (see Figure|2]l. 

The time averaged spectra of XISQ, XISl, XIS3, and 
HXD/PIN detectors are shown in Figure [3] The HXD/PIN spec- 
trum of 4U 1344-60 is background-dominated with the 27% of 

3 http://xmm.esa.int/sas/ 



the flux due to the source. This value is well above the back- 
ground systematic error (around 3-5%, see e.g. Suzaku-Memo- 
2008-03). The contamination of the 4U 1 344-60 HXD/PIN spec- 
trum by Centaurus B is negligible (see Sect. I7.2l in Appendix A). 
The 15-60keV flux of 4U 1344-60 is approximately 8 x 
10 _1 1 erg crrT 2 s _1 while the estimated flux from the XIS spectra 
of CenB is only 8 x 10~ 12 ergcirT 2 s _1 , i.e. one order of magni- 
tude lower. 

The underlying model in Figure[3]is an absorbe d power-law 
with the photon index T a 1.7. The tbabs model (I Wilms et al.l 
2000) was used for the photoelectric absorption. We fixed the 
column density of the absorption due to inter-ste llar material in 
our G alaxy to the val ue N H = 1.06 x 10 22 c m~ 2 ([Kalberlaet ail 
120051) . The survey by iDickev & Lockmanl (1 19901) gives a con- 
sistent value (Nh = 1.05 x 10 22 cnT 2 ). An additional absorption 
of comparable column density was required. We associated it to 
the host galaxy, although, regarding the very low Galactic lati- 
tude of 4U 1344-60, some fraction of it might be due to local 
variations in our Galaxy. The column density of the additional 
absorber was found to be Nh = 0.93 + 0.02 x 10 22 crrT 2 with the 
ZTBABS model. 

The only clear residuals from this simple model are around 
6.4 keV, and can be attributed to the iron fluorescent line. We 
added a Gaussian line to the model with the centroid energy 
and width as free parameters. The fit is then characterised by 
the C-value 7615 with 7312 degrees of freedom (v), or xl= 
5612/5452 » 1.03, respectively. The equivalent width of the 
line is 96+g 4 eV. The energy of the line was found to be £ = 
6.40+0.03 keV and its width cr = 0.085+°°^ keV (see Figure©. 

The corresponding FWHM velocity is « 8800 kms 1 , which 
suggests an origin in the outermost reg ions of an accretion disc 
or in a broad line region dPeterson 2006). T he optical m e asure- 
ments of a broad component of H a by Piconcelli et al.l ([2006) 
gives FWHM velocity twice smaller (» 4400kms~ I ) imply- 
ing that the iron line must originate closer to the black hole 
than the optical broad line region. However, compared to the 
XMM-Newton observa tion with E « 6.2 keV, and cr « 0.3 keV 
(IPiconcelli et al.ll2006l) the Suzaku iron line profile is much nar- 
rower. Therefore, we perform a simultaneous fit of both obser- 
vations in the next section. 
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Fig. 5. Comparison of Suzaku and XMM-Newton data residuals. Up: data residuals from the identical underlying model consisting 
of an absorbed power law and a Gaussian line. Bottom: data residuals from the final model containing an additional absorber with 
different covering fraction and slightly different slope of the power law for the XMM-Newton and the Suzaku observation. The front- 
illuminated spectra of XISQ and XIS 3 are jointly shown by black data points, the back-illuminated XIS 1 are red, and the EPIC-pn 
data are blue and emphasised by small circles. The data are binned for plotting purposes only. 



Table 1. Spectral analysis of 4U 1344-60 - the model parameters. 



model component 


parameter 


flat power-law model 


'part, cov.' model 


complex model 'fix' 


complex model 'free' 






XMM 


Suzaku 


XMM 


Suzaku 


XMM 


Suzaku 


XMM 


Suzaku 


local absorption 


«H,fully[10 22 cnT 2 ] 


0.86!°- 


0.93 ± 0.02 


0.92 ± 0.02 


0.95 ± 0.02 


0.98 ± 0.02 




«H,parl[10 22 CirT 2 ] 




7 5 +u 




8.2 ± 1.5 




cov. fraction [%] 






6±3 


50 ±3 


5±3 


41+15 


10 ±6 


power law 


r 


1 26 +004 

1 ' z ' u -0.02 


1-68^ 


1.66 ±0.09 


1.72 ±0.02 


i 79+O.O3 




1.83 ± 0.06 


reflection 


R 






0.48 ± 0.05 


2± 1 




-Efold 






100 (f) 


200 ± 100 




f [deg] 






45(f) 


85_ 6 




A^IO- 3 ] 


5 7 +03 


10.9 + 0.2 


7.1 ±0.4 


11.0± 0.3 


7 6 +0 - 7 

'•"-0.6 


11.2 ±0.2 


7.5 ±0.5 


12.2 ±0.6 




iVabstlO" 3 ] 




6±2 


7 +04 


7.6 ±0.3 


0.6 ± 0.4 


5±2 


1.4 ±0.8 


Gaussian line 


£ res , [keV] 


6 26 +o os 

"'^"-0.07 


6.40 ± 0.01 


6 29 +0 06 


6.40 ± 0.01 








cr [keV] 


ri 9^+0.09 
"••"-0.08 


0.08!™ 2 


91+0.06 

u ' il -0.07 


0.08 ± 0.02 








Mine [10- 5 ] 


10 + 3 


A 7+0-6 
^■'-0.4 


9±3 


4 7+0.5 

^■'-0.6 






fit goodness 


C/v 


9543/9110 


9484/9108 


9477/9114 


9460/9111 



Notes. Th e common parameter s for all the fits are the column density of the absor ber due to interstellar m atter in our Galaxy «h.ism = 1 -06 x 
10 22 cirr 2 iKal berla et alj|2005T) . and the cosmological redshift z = 0.012 ± 0.002 (Piconcelli et al. 2006). Solar abundances were assumed in 
the reflection models. The TVdi, parameter is a normalisation factor of the direct component while Mbs corresponds to the component affected by 
the partially covering absorber. They belong to the power-law, or the reflection model by PEXMON, respectively. The sign '(f)' means that the 
parameter was fixed to that value. 
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Fig. 6. Contour plot of the spin and the inclination angle mea- 
sured by the relativistic iron Ka line in the fiat power-law model. 
The individual lines correspond to 68, 90, and 99% confidence 
levels. 



4. Comparison with the XMM-Newton observation 

4.1. Evidence of long-term spectral variability 

Although the time-averaged observed flux did not change signif- 
icantly between the XMM-Newton and Suzaku observations the 
spectral shape is different. The simple power-law model with a 
Gaussian line applied to both data sets with the same parameters 
as from the Suzaku fit (besides the normalisation factors) does 
not provide an acceptable fit to the EPIC-pn spectrum with the 
total C/v = 10476/9112, or xl= 7340/6219 ~ 1.2, respectively. 
Figure [5] (upper panel) clearly shows that either the photon in- 
dex of the power law or the properties of the absorber changed 
between these two observations. 

4.2. Model with a variable power-law slope 

First, we allowed the column density of the local absorber 
to vary. The fit improved to C/v - 9877/9113, or xl~ 
6741/6220 ~ 1.08, respectively. The value of the column den- 
sity of the local absorber was found to be Nh,xmm = (1.5 ± 
0.1) X 10 22 cirT 2 . A significantly better fit (C/v = 9554/91 1 1, or 
xl— 6468/6218 « 1.04, respectively) was found when we also 
relaxed the power-law photon index. New values are Fxmm 
1 .24 + 0.04 and %mm = (0.85 + 0.06) x 10 22 cirT 2 . The fit 
goodness further improved by AC ~ 1 1 when, in addition, the 
iron line parameters were unbound. The resulting fit parameters 
are summarised in TableQ]in the column labelled "fiat power-law 
model". 

The energy of the iron line is significantly red-shifted to- 
wards i ts rest energy and broa dened, consistently with the find- 
ings bv lPiconcelli et al.l(l2006l) . We replac ed the Gaussian line b y 
the relativistic iron line model kyrline dDovciak et alj|20 04b). 
The resulting fit improved by AC « 22 with respect to the model 
with a broad Gaussian line, suggesting significant asymmetry of 
the line profile. The angular momentum of the black hole is pre- 
ferred to be high, i.e. the disc extends very close to the black 
hole. The inclination angle was found to be i = (22 + 7)°. We as- 
sumed the standard radial emissivity profile (i.e. reflection flux 
oc r~ 3 ) and isotropic angular emissivity. Figure [6] shows the con- 
fidence levels of the spin and the inclination angle. 



Although this model can satisfactorily explain the XMM- 
Newton spectrum, the narrow line from the Suzaku spectrum 
suggests that the reflection originates much further from the cen- 
tral black hole. A major part of the iron-line profile can be indeed 
explained by the same narrow line as seen in the Suzaku spec- 
trum (see Sect. l5.2l for more discussion). Therefore, we proposed 
an alternative model explaining the spectral change between the 
XMM-Newton and Suzaku observation. 

4.3. Model with a complex absorber's variability 

iPiconcelli et alj d2006l) suggested that a partially covering ab- 
sorption might be responsible for measuring the fiat spectral 
shape in the XMM-Newton observation. Employing the model 
with more complex absorption did not improve the statistical 
goodness of the fit, but yielded a photon index m ore similar to 
the v alues measured in Seyfert galaxies (see e.g. iBianchi et al.l 
2009). Therefore, we also included the partially covering absorp- 
tion in our model, first only to the XMM-Newton data. We bound 
the iron-line parameters between the Suzaku and XMM-Newton 
spectra for the case that the spectral curvature at around 5-6 keV 
could be explained by partially absorbing cl ouds, similar l y as in 
sp ectroscopy st udies of different sources bv lMilleret all (2008) 
or lTurneretai] (120111) . 

The photon index of the XMM-Newton spectrum, F = 1.6 + 
0.1, is indeed found to be more consistent with the Suzaku obser- 
vation. The goodness of the fit is also better than in the previous 
case with C/v = 9499/9109, or^ 2 = 6408/6216 * 1.03, respec- 
tively. The column density of the fully covering absorption was 
found to be A^xMM.ftdiy = 0.86 ±0.11 x 10 22 cm~ 2 . The fraction 
of the partially absorbed power law is 44 + 12 % with the column 
density Afc,xMM, P art = 7.1+ 2 ;° x 10 22 cirT 2 . 

As a next step, we included the partially covering absorp- 
tion also for the Suzaku data. We simply assumed that only the 
fraction of the partially covering absorption changed between 
the times of the two observations (allowing the column den- 
sities to vary did not improve the fit significantly). The col- 
umn density of the fully covering absorption was found to be 
Awfully = 0.92 + 0.02 x 10 22 crrT 2 , the column density of the 
partial covering A^/, P art = 7. 2 x 10 22 crcT 2 , and the frac- 
tion of the partially absorbed power law 44*\l % for the XMM- 
Newton and 5 + 3 % for the Suzaku spectrum. The fit goodness 
is C/v = 9490/9110 ,or xl = 6347/6217 a 1.02, respectively. 

IPiconcelli et al.l (120061) rejected the model with the partially 
covering absorber (their "model E"), because it did not explain 
all the residuals at the iron-line energy band. Indeed, when we 
unbound the iron-line parameters, we obtained similar results 
as with the fiat-powerlaw model. The energy of the line is found 
red-shifted, E = 6.3+0.1keV, and broadened, cr = 0.2 + O.lkeV. 
The equivalent width of the line is £Wsuz = 96^2 eV with 
Suzaku, and EW X mm = 160^ eV with XMM-Newton. The dif- 
ference consists in the width of the line. A more detailed look on 
the iron-line residuals is shown in Fig. [7] The partially-covering 
absorber does not explain why the iron line is broader and more 
redshifted in the XMM-Newton spectrum. However, we consider 
this model further, as it is more consistent with the Suzaku mea- 
surements than the flat power-law model. 

4.4. More complex reflection models 

The iron line detected in the XIS spectra implies the presence of 
the cold reflectio n component in the spectrum. We used further 
PEXMON model dNandra et all l2007h . which combines the re- 



5 



J. Svoboda et al.: Active galaxy 4U 1344-60: did the relativistic line disappear? 



o 



o 
E 




Energy [keV] 

Fig. 7. Spectral residuals in the iron line energy range. The resid- 
uals are calculated against model consisting of power law with 
the partially covering absorber. Black data points correspond to 
the spectra of front-illuminated Suzaku/XIS detectors. The red 
ones (with small circles) correspond to EPIC-pn spectrum. The 
data are binned for plotting purposes only. 
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Fig. 8. Final model consisting of a partially covering absorber 
and cold reflection. The model is shown by red or black re- 
lated to the XMM-Newton or the Suzaku observation, respec- 
tively. The components affected by the partially covering ab- 
sorber are shown by dashed lines while the direct components 
by solid lines. 



flectio n continuum by PEXRAV model ( Magdziarz & Zdziarskil 
1 1995b with the appropriate strength of the fluorescent lines of 
iron and nickel. Because the iron line detected by XIS detec- 
tors is narrow, the reflection must occur sufficiently far from the 
black hole and is rather unlikely to be variable between the two 
observations. Therefore, we simply assumed the same reflection 
model for both Suzaku and archival XMM-Newton data. 

Neither the inclinati on angle nor the origin of the reflection 
is known in this galaxy. iPiconcelli et alJ 12006) determined the 
type of Seyfert 1.5 galaxy from the optic al measuremen ts. In 
terms of the AGN unification scenario bv lAntonuccil d 19931) . this 
would correspond to an intermediate value of the inclination, 
i.e. i w 45 degrees. Therefore we considered two cases: first one 
with keeping the values of the inclination and the folding energy 
fixed in the fit (called as "complex fix" in Table [T]) and second 
one with these parameters free ("complex free" in Table [1] or 
"final" model hereafter). 

The fit goodness of the final model is characterised by C/v — 
9460/9111, or xl= 6316/6218 ~ 1.02, respectively. The final 
model is plotted in Figure [8] in the 1 - lOkeV energy range. It 
is demonstrated there how the fraction affected by the partially 
covering absorber decreased between the two observations. The 
iron and nickel fluorescent lines as well as the iron edge are the 
properties of non-relativistic cold reflection. 

Although the spectrum is described rather satisfactorily by 
the final model, the residuals at the energy range around 6 keV 
still persist in the XMM-Newton spectrum. The data/model ratio 
is shown in the bottom panel of Fig. [5] In order to demonstrate 
the significance of the remaining residuals, we added a Gaussian 
line to the final fit. The fit improves by AC = 22 (or Ax = 23, re- 
spectively). This is a statistically significant result as the random 
improvement should be of order of AC ~ 6 for the three extra 
parameters added to the model. The energy of the Gaussian line 
was found to be E = 6.15+°} 2 keV, the width cr = 0.21+°J 2 keV, 
and normalisation factor A^jauss = 5.8^ \ X 10~ 5 . The equivalent 
width of this feature is 90*12° e V. This additional line can be at- 
tributed to a redshifted iron line emitted from inner parts of the 
accretion disc - see Sect. l5.2l for more discussion on this spectral 
feature. 



5. Discussion 

5. 1 . Presence of a partially covering absorber 

The Suzaku observation of 4U 1344-60 reveals that the spectrum 
can be characterised by a photon index of the p ower law T » 
1.7, w hich is a typical value for Seyfert galaxies (Bianchi et al. 
2009). The different slope of the XMM-Newton spectrum can 
be explained by a very low photon index (F < 1 .3) or by the 
presence of a partially covering absorber. A statistically better 
fit was obtained with the model with the partially covering ab- 
sorber. Moreover, this model yields flux at high energies more 
consistent with the measured values by the Integral and Swift 
satellite, as demonstrated in Table [2] 

High-energy flux predicted by different models (flat power 
law, model with the partially covering absorber, final model) is 
compared with the values measured by the Integral and Swift 
satellitfl Better agreement is obtained for the scenario with the 
partially covering absorber (best for the final model). However, 
this cannot be regarded as a final proof owing to a possible 
long-term variability between the non-simultaneous observa- 
tions. Moreover, a low high-energy cut-off in the flat power law 
would significantly decrease the predicted hard X-ray flux by 
this model. 

An additional fully-covering absorber with Nh ~ 10 22 cm~ 2 
along our line of sight to the nucleus of 4U 1344-60 was required 
by the data. It may be identified with dust lanes or star-forming 
regions in the i nnermost h ost galaxy disc and not associated with 
the AGN torus (Matt 2000). A very larg e contamination from the 
host galaxy has also been claimed by IVasudevan et al.l (1201 Ol) 
in order to explain the infrared spectral energy distribution of 
4U 1344-60. This scenari o is consistent with the i ntermediate 
type classification (see e.g. lMaiolino & Riekdll995|). which w as 
originally proposed for 4U 1344-60 by IPiconcelli et a l. (2006). 



4 A power-law model was used to measure Swift- 
BAT flux from the 58-month hard X-ray survey data. 
http://heasarc.gsfc.nasa.gov/docs/swift/results/bs58mon/ 
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Table 2. Flux of 4U 1344-60 measured by different instruments and predicted values by various models. 



instrument (- model) 


flux [10" n erg cm" 2 s" 1 ] 




2 - lOkeV (observed) 


2 - lOkeV (intrinsic) 


20-40keV 


40- lOOkeV 


20- lOOkeV 


XMM - flat T 


4.16 + 0.04 


4.71+0.04 


7.7 + 0.1 (p) 


19.0 ± 0.2 (p) 


26.7 ± 0.3 (p) 


XMM - part. cov. 


4.04 ± 0.04 


5.56 ± 0.04 


4.3 + 0.1 (p) 


7.5 + 0.1 (p) 


1 1.9 + 0.2 (p) 


SUZAKU - part. cov. 


3.83+0.01 


4.58 ± 0.02 


3.29 ± 0.02 


5.49 ± 0.03 (p) 


8.75 ± 0.05 (p) 


SUZAKU - final 


3.83+0.01 


4.67 ± 0.02 


3.35+0.02 


4.09 ± 0.03 (p) 


7.44 ± 0.05 (p) 


Integral 






3.2 + 0.1" 


4.1+0.2° 


6.6* 


Swift-BAT 






2.7 + 0.1 


4.0 + 0.2 


6.7 + 0.3 



Notes. A sign (p) means that the value was predicted by the model from lower X-ray energies. (a> Bird et al. (2010), (6) Pan essa et al] 12008). 



5.2. The iron line in the XMM-Newton observation 

Model consisting of a power-law component modified by the 
fully and partially covering absorber and the narrow iron fluores- 
cent line represents a good simultaneous fit to both, the Suzaku 
and the archival XMM-Newton spectrum, apart from the emis- 
sion residuals at the moderately redshifted iron line energy band 
that are entirely related to the XMM-Newton spectrum. Figure [7] 
shows this energy band for the Suzaku-XI S and XMM-Newton 
spectrum in detail. Piconcelli et al] ([2006) concluded from the 
XMM-Newton spectrum that the total iron line profile is most 
likely formed in the innermost accretion disc. However, the line 
profile measured by Suzaku is much narrower. Therefore, we in- 
terpreted the XMM-Newton iron line profile as the narrow com- 
ponent (the same as in the Suzaku spectrum) plus an additional 
enhanced emission from the innermost accretion disc. 

Such an enhanced emission can be explained, e.g., by a 
temporary irradiation of a part of the accretion disc by a flare, 
for instance due to a magnetic reconnection in the corona 
dCzernv et ai1l2004t lUzdenskv & Goodmanll2008l) . With the en- 
ergy E ss 6.1keV, it is significantly redshifted when asso- 
ciated to the iron fluorescent line. Therefore, the probable 
origin is from the innermost accretion disc where the radi- 
ation is affected by strong gravitational redshift. A transient 
emission feature at a simila r energy was report e d for another 
Seyfert galaxy NGC3516 dTurner et al] l2002t iBianchi et all 
2004). Much longer exposure time allowed to trace the temporal 
evolution of the emission and reveal possib le periodicity d ue to 
an orbiting spot in the accretion disc dlwasawa et al]|2004l) . 

We performed several tests using KY model dDovciak et al] 
2004b) including also non-axisymmetric models (due to the 
presence of a partially covering absorber). We tried to split the 
XMM-Newton observation into several parts and look for the 
change of residuals. We also generated a light curve related only 
to the redshifted component of the iron line and compared it 
with the total light curve. However, the insufficient data qual- 
ity of the short XMM-Newton observation does not allow us to 
uniquely constrain neither the geometrical structure of the emis- 
sion region nor its temporal behaviour. It was not possible to dis- 
tinguish between a "ring-shaped" static emitting region and an 
orbiting spot. The evidence of such spectral features, however, 
satisfy the studies o f a locally enhanced emission from parts o f 
the accretion disc (Dovciak et al. 2004a; Goosma nn et al l2007h . 



and will be potentially r esolved with future h igher resolution X- 
ray instruments (see e.g. lSochora et al]|201 lb . 

The absence of any relativistic features in the Suzaku spec- 
trum is rather puzzling. Although the 4U 1344-60 observa- 
tion revealed a complex absorption the Suzaku observation pro- 
vided us the clearest view on the nucleus. A possible explana- 
tion for non-detection of the redshifted iron-line emission could 
be due to a truncation of the disc at a further radius. A trun- 
cat ed accretion disc was reported for several other sources (see 
e.g.lMatt et al]l2005HMarkowitz & Reevesll2009HSvoboda et all 
2010). The physical reason for such a truncation might be due 
to a very low accretion rate when the accretion flow would be 
advection domi nated and would not form a thin accretion disc 
(Esin et alii 19971) . or due to an over-ionisatio n of the innermost 
part of the accretion disc by strong irradiation dDone et all 2000; 
iBhavani & Nandrall201 ll) . 

The change of the iron line from a relativistically broadened 
profile to a narrow one as a likely result of the stat e transition 
of the source, was reported for a quasar Q0056-363 (Matt et al. 
2005). Different photon indices were measured between the two 
observations but on contrary to our case the flatter index was 
related to the narrow line detection. An explanation in terms of 
the disc truncation is therefore unlikely in the case of 4U 1344- 
60 since it does not suite well in a truncat ed disc scenario in the 
hard states of Galactic bla ck hole binaries (Lubihski & Zdziarskil 
120011: lFenderetaljl2004h . Moreover, the measured flux during 
the Suzaku observation is very close to that measured by XMM- 
Newton. The accretion rate de rived from the infrared and X-ray 
observations is X » 0.03/lEdd dVasudevan et alJ[2010h . which is 
somewhat typical value for other Seyfert galaxies as well. 

Therefore, the ionisation may represent a more plausible ex- 
planation. The illuminating flux is a strongly radially dependent 
functi on, especially if the coron a is very comp act (Matt et al. 
[19931: iRevnolds & Fabianll2008t ISvoboda et all 120121) . The in- 
nermost disc can be over-ionised owing to the combination of 
sufficiently high irradiation flux and relatively low density of the 
accreting material. The ionisation state could, however, locally 
decrease if the density is locally increased (e.g. due to a shock 
wave), or if some clouds partly obscure the illuminating source. 
The latter one could happen during the XMM-Newton observa- 
tion. As a result, some additional fluorescent emission was de- 
tected from the innermost part of the accretion disc. 
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Fig. 9. Contour plot between the inclination angle and reflection 
parameter of the final model. The individual contours correspond 
to 68, 90, and 99% confidence levels. 



An alternative explanation for the missing imprints of the in- 
nermost region could be due to a distant corona. Speculatively, 
if the scattering medium is associated with the partially cov- 
ering absorber this would also explain the change between the 
XMM-Newton and the Suzaku observation. The decreased cov- 
ering fraction would be due to an out-flowing corona and thus, 
less emission from the innermost part of the disc would be de- 
tected in the Suzaku observation. Smaller intrinsic X-ray flux 
(see Table [2]l during the Suzaku observation would match well 
this scenario. 

5.3. Reflection models and the angle of inclination 

The reflection model PEXMON was used for a consistent descrip- 
tion of the iron Ka, the iron edge and the Compton hump. The 
reflection fraction R, the folding energy Ef n and also the in- 
clination angle i are, however, not constrained very well and 
suffer from a mutual degeneration of the parameters. Moreover, 
the link between the Compton hump and fluorescent lines might 
not be necessarily fulfilled. Especially in the case of reflection 
from the torus a different geometry should be taken into account 
(iMurphv & Yaqoobl l2009t lYaqoobl l201~2h . We also performed 
several test fits with the PEXRAV model and a Gaussian line sep- 
arately and we get almost identical results like with PEXMON, 
i.e. favouring a high inclination angle measured relatively to the 
normal of the orbital plane. 

However, it is necessary to stress that the inclination angle 
intervenes only as a projection effect in the PEXRAV or PEXMON 
model. The shape of the emission is independent of the emission 
angle and the total spectrum is averaged over all angles. As a re- 
sult, both R and i affect only the normalisation and are therefore 
degenerated. This is apparent in Figure [9] where the confidence 
levels of the reflection parameter R and the inclination angle i of 
the final model are shown. Illumination of the accretion disc by 
an isotropic source on top of it corresponds to the value R — 1 . 
The larger values of R would suggest that the more observed 
radiation comes from reflection than from the direct component. 

The resulting high inclination angle of the final model, which 
is also suggested for the case R — 1, is likely a consequence of 
the low contrast between the Compton scattering hump and the 
primary continuum, the former being determined by the combi- 
nation of strong Compton scattering and low photoelectric ab- 



sorption cross-section. The lower inclination would imply the 
Compton hump to be more prominent given the observed equiv- 
alent width of the iron line. 

This would would imply two possible scenarios. First, the 
system could be oriented almost face-on but we see the reflection 
on the distant broad line region or the torus, which gives a com- 
plementary reflection angle. This is s upported by the detect ion of 
broad lines in the optical spectrum dPiconcelli et al.l l2006). It is 
also consistent with low inclination angle of the innermost accre- 
tion disc measured by the relativistic iron line model in the flat 
power-law scenario for the XMM-Newton spectrum (see Fig. [6j>. 
The second scenario assumes that the reflection occurs at the 
more distant parts of an accretion disc which is oriented almost 
edge-on relatively to the observer. 

The edge-on disc hypothesis is not in a strong contradiction 
with the optical measurements. The standard unification sce- 
nario has been ques tioned in several sources, for a review see 
iBianchi et alJ (120121) and references therein. The dumpiness of 
the absorber is supposed to be another important parameter in 
the proper classification and the type of AGN could not be en- 
tirely determined by the inclination angle but also by a prob- 
abilit y of the absorption clouds i ntercepting the line of sight 
(Nen kova^alfcOOalEih^l^OllI) . This is supported by a large 
variability of the a bsorber observed in several sources (see e.g. 
iRisaliti et aTll2005l) . 4U 1344-60 has also revealed partially cov- 
ering absorber with a significant change between the archival 
XMM-Newton and recent Suzaku observation, and can be thus 
another example advocating for an updated unification model. 

Moreover, the inclination angle of the acc retion disc and 
the h ost galaxy might be misaligned (see e.g. [Guainazzi et al.1 
1201 ll) . An evident misalignment of the circum-nuclear mat- 
ter with the radio jet, which is naturally b elieved to li e in th e 
black hole rotation axis, was reported by Rab an et all {2009). 
IRisaliti et al.l (1201 ll) analysed the distribution of the [O m] emis- 
sion line equivalent width in order to constrain the inclination in 
a large sample of AGNs, and reported several sources with edge- 
on discs but not heavily absorbed by a circum-nuclear matter. 4U 
1344-60 could belong to this kind of sources as well. 



6. Conclusions 

The X-ray continuum spectrum of the Seyfert galaxy 4U 1344- 
60 is dominated by a power law with a standard value of the 
photon index 1 .7 +0. 1 . A narrow fluorescent line with equiv- 
alent width « 100 eV is an apparent feature in the XIS spectra. 
The width of the iron line is cr = (0.08 ± 0.02) keV, i.e. * 9000 
km/s of FWHM, suggesting its origin from a distant matter not 
affected by relativistic smearing. The outermost parts of an ac- 
cretion disc, broad line region clouds or torus are possible can- 
didates for the reflector. Its profile is, however, significantly nar- 
rower than in the archival XMM-Newton observation. The de- 
tected red-shifted iron line emission during the XMM-Newton 
observation was likely a temporary feature either due to locally 
enhanced emission or decreased ionisation in the innermost ac- 
cretion flow. 

The spectral shape has also changed between the XMM- 
Newton and Suzaku observation. We interpret this change as due 
to the presence of a partially covering absorber. The absorbed 
fraction of the primary radiation has decreased from around 45% 
to less than 10%. The resulting model is then characterised by 
the same power law slope. It is also consistent with the flux mea- 
sured (non-simultaneously) by the Integral and Swift satellite at 
high X-ray energies. The spectral variability makes 4U 1344-60 
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an interesting target for further examination within a monitoring 
programme by an X-ray satellite. 
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7. Appendix A 

7. 1 . Timing and spectral properties of Centaurus B from 
ASCA to Suzaku 

Centaurus B is one of the nearest and brightest radio galaxies, 
at cosmologica l redshift z = .01215 and flux density « 250 Jy 
at 408 MHz dJones et al .11200 ll) . It is located only at 13' from 
4U 1 344-60 on the sky. The proximity of the two sources did 
not allow us to set an observing configuration in which we can 
exclude Centaurus B from the HXD detectors. Therefore, we re- 
quired an accompanying lOks observation targeted directly on 
Cen B to estimate its X-ray flux in HXD/PIN energy range. Here, 
we present our spectral results from this short Suzaku obser- 
vation using the XIS detectors together with the analysis of its 
long-term variability. The data reduction procedure is the same 
as described in Section|2] 

Historically, the flux of Centaurus B was found to be only a 
factor of a few smaller than 4U 1344-60. The observed flux in 
the 2-10 keV was 7.3 + 0.2 x 10" 12 er g cm" 2 s" 1 in the archival 
ASCA observation dTashiro et alJI 19981) . and 5.6 ± 2xl0~ 12 erg 
cm" 2 s" 1 in the XMM-Newton observation (see also lTashiro et al.l 
120051) . i.e. around 6-7 times fainter than 4U 1344-60. The 
Integral satellite measured hard X-ray fluxes of /20-40keV = 
(0.76 ± 0.08) x 10" 1 1 erg cm" 2 s" 1 an d Ao-iookev = (1.1 ± 0.2) x 
10 erg cm V 1 (iBird et alJl2010h . around 4 times less than 
the 4U 1344-60 fluxes. 

Figure [10] shows X-ray light curves from the archival ASCA 
and XMM-Newton observations. The light curves look stationary 
with no trend in the flux. We fitted the light curves by polynomial 
fits, and we found that a constant function represents the best 
fit. Most of the apparent variability in the light curves is due to 
Poisson noise whose level is shown by the dashed lines in the 
plots. These results clearly suggest very low variability and, in 
particular, no trend in the average flux of CenB on a w lOOks 
time-scale. 

The 2-10 keV flux measured with the Suzaku/XIS detec- 
tors is (3.9 + 0.2) x 10" 12 erg cm" 2 s _1 , which follows a con- 
tinuous decrease of the flux since the ASCA observation (see 
Table|3]l. The extrapolated 15-60 keV flux is then approximately 
8 x 10" 12 erg cm" 2 s" 1 , which is one order of magnitude less than 
4U 1 344-60. Thus, in the hard X-ray regime, the flux of 4U 1 344- 
60 clearly dominates. 

Although the X-ray brightness of CenB has decreased by 
a factor of two from the archival ASCA observation in 1995 to 
the recent Suzaku observation in 2011 the spectral shape has 
not changed (see Table [3} and can be characterised by a sim- 
ple power law with the photon index « 1.6 and local absorption 
with the column density Nh ~ 0.7 x 10 22 cm" 2 . The Suzaku/XIS 
spectrum of Centaurus B is shown in Fig. [TTJ The underlying 
model is TBABS*ZTBABS*ZPOWERLW in the XSPEC notation 
with the best-fit values from the Table [3] The statistical good- 
ness of the fit can be characterised by C/v = 177/162 (or 
X 2 /v = 169/162 « 1.04). Some residuals from the model are 
apparent in the spectrum but they are likely due to a less than 
perfect subtraction of the background, considering the relatively 
low flux of CenB (the fraction of the background flux is around 
30%). 



9 



J. Svoboda et al.: Active galaxy 4U 1344-60: did the relativistic line disappear? 



Z5 

o 
o 




£L 0.8 







CO 



o 
o 



0.6 







50 



25 



30 



100 150 200 5 10 15 20 

time [ks] time [ks] 

Fig. 10. Light curve of Centaurus B. Left: ASCA (200 ks). Right: XMM-Newton (30ks). Horizontal solid lines represent the average 
flux. Dashed lines show the 3cr levels where <x is the expected standard deviation due to Poisson noise. 



Table 3. Evolution of the X-ray flux and spectral slope of 
Centaurus B. 





ASCA 


XMM-Newton 


SUZAKU 


2- lOkeV flux 
[10~' 2 ergcnr 2 s -1 ] 


7.3 ± 0.2" 


5.6 ± 0.2 


3.9 ± 0.2 


photon index 


1.64±0.07 fl 


1.6 + 0.1 


1.6 + 0.2 


column density* 


1.76±0.11 fl 


1.7 + 0.1 


1.7 + 0.3 



(a) Tashiro et al., 98 
{b) including Galactic absorption (N H = 1.06 x 10 22 cirT 2 ) 



7.2. Estimate of the contamination of the 4U 1344-60 
HXD/PIN spectrum by Centaurus B 

Although the high energy flux of Centaurus B is one order 
of magnitude lower, we considered it for the proper cross- 
normalisation between XIS and HXD/PIN instruments. We per- 
formed a fit of HXD/PIN spectrum of 4U 1344-60 with a sin- 
gle and double power-law model. The second power-law com- 
ponent did not improve the fit at all. However, a slightly flatter 
photon index was preferred using the single power-law model. 
This would be an expectable influence of the contamination 
by a harder Centaurus B spectrum with the photon index F « 
1.6 + 0.2. Although the determined error of the photon index 
from the short Suzaku observation is rather high, the value of 
1 .6 has been consistently measured since the earliest ASCA ob- 
servation (see Table |3). 

We used the best fit values obtained from the analysis of the 
XIS spectra of Centaurus B multiplied by a constant between 
and 1. We fitted the XIS and HXD/PIN spectra simultaneously 
with the final C-value 7847 with 7505 degrees of freedom, or 
x\- 5839/5645, respectively. This means only a marginal im- 
provement with AC « 3 (or Ax ~ 4, respectively) compared 
to the model without any contamination taken into account. The 
best-fit value for the constant is c — 0.1^2"*. Given the fact that 
the flux of Centaurus B is one order of magnitude weaker than 
for 4U 1344-60 the contamination is 5% at maximum. 
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Fig. 11. X-ray spectrum of Centaurus B observed with Suzaku 
XIS0 (black), XIS 1 (red), XIS 3 (green). 
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